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Abstract 

We study quantum dynamics of three dressed Bose-Einstein condensates in a 
high-Q cavity. The quasiparticle excitation spectrum of this system is found 
numerically. The stability of the quasiparticle excitation is analyzed. It is 
shown that there exist instabilities in the excitation spectrum. 
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Recently, multicomponent Bose-Einstein condensates (MBEC's) have attracted 

much attention due to the recent experimental realizations of two and three component 
Bose-Einstein condensates (BEC's) [§],|J. MBEC's offer new degrees of freedom, which give 
rise to a rich set of new phenomena that do not exist in a single condensate. Especially, 
the interaction of BEC's with light leads to fascinating effects including the extreme slowing 
down of the speed of light || and matter- wave four- wave mixing ||. Goldstein, Wright and 
Meystre (GWM) [0] proposed an approach to produce MBEC's, called dressed BEC's, inside 
a high-Q multimode cavity, and investigated the quasiparticle excitation spectrum for the 
case of two dressed BEC's which dress one-photon. The purpose of this letter is to study 
the quasiparticle excitation spectrum for a system of three dressed BEC's which dress two 
photons. We will investigate quantum dynamics of the three dressed BEC's, and calculate 
numerically their elementary excitations. 

Consider the GWM model of dressed BEC's [7], which consists of a condensate which 
interacts with two counter-propagating modes by a high-Q ring cavity with n\ photons 
being in the mode 1 and n 2 photons being in the mode 2. We consider a state com- 
posing n = fix + n 2 photons and N atoms which comprises the condensate confined by 
an external trapping potential U(r). The state of the system can be written \^N,n) = 
E ni ,n 2 / dr i-- I dr N f nuna (r x , r N , t)^{v l )...^{r N ) 

x |0, rii, n 2 ), where ^(r,) are the atomic field operators, / ni) n 2 ( r i> r N, t) is the many- 
particle Schrodinger wave function for the condensate, which can be factorized as a 
product of Hartree wave functions 4> ni ,n 2 representing states which the atoms occupy, 
fm,n 2 (ri,...,r N ,t) = l\f=i 0m,n 2 (ri, t). 

Let us consider such a trapping potential || U(r) = mur[(r5_ + X 2 z) 2 , where r = (r±, z), 
rj_ being the transverse position coordinate, u± the transverse angular frequency of the trap, 
and A = uj z /u>l is the ratio of the longitudinal to transverse frequencies. For the cigar struc- 
ture condensate ||, it is reasonable to assume that A <^ 1, and the transverse structure of 
the condensate is determined by the ground state solution of the transverse potential i>(rj_). 
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Then the Hartree wave function can be expressed as ni! „ 2 (r, t) = v(r±)e~ lU} - Lt ip nijn2 (z : t). 
From the Hartree variational pronciple and after integrating over the transverse coordinate 
r_i_, one obtains Gross-Pitaeviskii equations (GPE's) 

#ni,na(f,T) = H L l/j nun2 (£, t) 

+V\ ^ni,n 2 (C,r) \ 2 ^ ni ,n 2 (^r) 

+g[\jni{n 2 + l)^ ni _i >na+ i(^, r) 

(C,r)], (1) 

where = —d 2 /dt; 2 + £ 2 /4, with £ = z/a z , a z = ^Jh/2mu z , r = u z t, g and r\ are the linear 
and nonlinear coupling coefficients, respectively. 

The dressed BEC's are the eigenstates of Eq.(l), and are quantum superposition states of 
states with different photon numbers (ni,n 2 ). We consider the case of three dressed BEC's. 
In this case, only the states with (2,0), (1,1) and (0,2) are concerned. The coupled GPE's 
(1) reduce to 

#2,0 (f, r ) = H L ip 2 ,o(^ r ) + v / 2#i,i(^ t) 

+ V\ ^2,o(£,T) |V2,o(£,T), 

#i,i(£,t) = H L rl> Xil (£,T) + y/2g[rl) 2i0 <£,T) 

+ +ti\M£,t) | 2 ^i,i(e,r), 

#o,2(^,r) = i/ L ^o,2(^r) + v / 2#i,i(^,r) 

+ v \^ 2 (Z,t) |Vo,2(e,r). (2) 

The exact analytic stationary solutions are generally not available for the coupling GPE's 
of the dressed BEC's. We here consider their solutions under the Thomas-Fermi approxi- 
mation (TFA), in which the non-linear interaction term dominates over the kinetic-energy 
term [9]. Setting i/) ni! n 2 (£, t) = e~ 4/1T #ni,n 2 (0 with /i being the chemical potential, under the 
TFA, coupling GPE's (2) reduce to 

/Ao(0 = |fM0 + V2g9 1A (0 + rj\ 2>o (fl | 2 2)O (O, 
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/*M0 = ^MO + V2g0 ltl (t) + v I do,2(0 I 2 M£), 

/A,i(0 = |£ 2 M0 + ^#0,2(0 + 02,0(0] 

+>HM0l 2 M0- (3) 

For simplicity, we consider only solutions of two cases for 6*2,0 (0 an d 6*0,2 (0 : ( a ) Out-of- 
phase solution with #2,0 (0 = — #0,2 (0; O 3 ) In-phase solution with 62,0 (0 = #0,2 (0- 

(a) Out-of-phase solution. In this case, #2,0 (0 — — $0,2 (0, from Eqs.(3) we obtain the 
following dressed state nonzero solutions 



MO^/t/. M0 = «, 4/i>e 2 . (4) 



Let £ m be the boundary of the dressed condensates at which the TEA solutions vanish. 
From Eq.(4) it follows that the chemical potential /i = The value of £ m is determined by 
the normalization condition of the condensate wave function J2m,n 2 f-^ m ^£ I ,712(0 | 2= 1, 
which leads to £ m = [3r]/2] 1 ^ 3 . So we get the following out-of-phase solution 

I MO l 2 H MO l 2 = ^jp, I MO l 2 = 0. (5) 

It is interesting to note that the out-of-phase solution of the three dressed BEC's are the 
same as those of two dressed BEC's found in Ref. [0] due to the vanishness of 6*1,1(0- This 
reflects the fact that the three dressed BEC system can transit to a two dressed BEC system. 
Physically, this is reasonable, because a three dressed BEC system naturally becomes a two 
dressed BEC system when one of the three dressed condensates vanishes. 

(b) In-phase solution. In the case, 6*2,0(0 = $0,2 (0; Eqs.(2) reduce to two equations: 

(n - ~f)M0 = ^01,1(0 + v I MO I 2 MO. 

+v 1 oi,i (o r mo, ( 6 ) 

which can be solved numerically. Assume £ m to be the boundary of the dressed condensates, 
from Eqs.(6) we then obtain the chemical potential ^ = ±2g + £^/4. Here the value of £ m 



depends on the normalization condition J2 ni ,n 2 f-£ m d£, | ^ni,n 2 (0 | 2= 1- For example, when 
g = 7.29 and rj = 80, we numerically find that = 4.584 and £~ = 4.430 coresponding to 
the chemical potential being /i + and ji~ , respectively. 

We now find the elementary excitations of the system by linearizing Eqs.(6) around 
the dressed state solutions: ?/> rai>ri2 (£, r) = e _i ' tT [0 ni>n2 (O + Uni,n 2 (0 e ~ <h,T + v* un2 (Z,)e iuJT }, 
where(ni,n 2 ) = (2,0), (1, 1) and (0,2), u(£) and v(£) represent small perturbations arround 
the dressed state with energies /i±ui. Substituting the expression of ip ni ,n 2 (^ T ) m t° Eqs.(6), 
in first order, we get the system of six equations for the linearized perturbations u(£) and 
v(0: 

um 2t0 (0 = [H L + 2 V | 6 2fi (0 | 2 -/xKo(0 

+V2^Ui,i(0+^2,o(OVo(0, 
uu 1 , 1 (Z) = [H L + 2ri\ 9 hl (0 | 2 -/i] Mlil (0 

+V2g[u 2 , (O + «o, 2 (0] + ^i(0Vi(0. 

uuo i2 (£) = [H L + 2 V I O , 2 (O | 2 -n]uo, 2 (.0 
+y/2g Ul>1 (0 +77^0,2(0X2(0, 
-uw 2 ,o(0 = [if L + 2 V I 2)O (O I 2 -*,o(0 

+^l,l(0+^2,o(OVo(0> 
= [#L + 2r/| 1;1 (O | 2 - /i]^l,l(0 

+^h, (0 + ^0,2(0] + ^l,l(0V,l(0, 

-uwo, 2 (0 = [H L + 2t7| o , 2 (O I 2 - ax]u 0) 2(0 

+v^i,i(0+^o, 2 (OV), 2 (0- (7) 

The normal modes of the above coupled equations are identical to the elemantary exci- 
tations determined via Bogoliubov method. In order to obtain the elementary excitations, 
we expand them in terms of eigenfunctions of the trapping potential ^(0 as 

*V,n 2 (0 = ZXi,n 2 ft'(0> 

v 
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Uf) = EC^)- ( g ) 

V 

where $,(0 satisfies H L q v ({) — (u + §)g„(0- 

It is difficult to get an exact solution for the normal modes of Eqs.(7). In what follows, 
we give numerical solution of them by using the TFA. For the out-of-phase solution (4), we 
note that | #2,0 (0 f — I #0,2 (£) | 2 — I #2,o(£ = 0) | 2 ^ for normal modes localized close 

to the center of the trapping potential due to £ m ^> 1. 

Substitution of Eqs.(8) into Eqs.(7) gives a system of linear equations for the coefficients 
K,i,n 2 an< ^ c ni,n 2 - The nonzero solution condition of the system of linear equations gives rise 
to the excitation spectrum of the elementary excitations. In the Fig. 1(a) we have plotted 
the excitation spectrum in terms of the index v of the linear oscillator mode for the out-of- 
phase solution when 77 = 80, g = 7.29 and rj m = 4.9324. Fig. 1(b) is the enlargement of the 
part of —10 < uo 2 < 20 for the curve c. From Fig. 1(a) we see that in the region of < v < 2, 
there exists one type of excitation spectrum, which is stable due to uj 2 > 0. In the region of 
v > 2 we can observe three types of excitation spectrum. Two of them indicated by curve 
a and b are stable excitations. The third one indicated by curve c is instable in the region 
of 14 < v < 18 since u 2 is negative. This is explicitly indicated in Fig. 1(b) through the 
enlargement of the curve c. 

Similarly, we can numerically express the excitation spectrum for the in-phase solution 
in terms of the index v of the linear oscillator mode. In Fig. 2, we give rise to the excitation 
spectrum for the in-phase solution when r\ = 80, g = 7.29 and i] m = 4.932, where we have set 
that I 2 ,o(O I' = I MO I' - I M£ = 0) | 2 w 0.0336, | | 2 ~ | 9 hl (£ = 0) | 2 w 0.0925 

for the branch and | 2jO (O | 2 =| MO l 2 ~ 0.04925, | M (O | 2 = | M (O) | 2 w 0.0775 
for the branch which lead to = 4.584 and £~ = 4.431, respectively. From Fig. 2(a), 
we see that there exist different excitation spectrum for three different regions of v. In the 
first region < v < 8, we can find three kinds of excitation spectrum, and all of them are 
stable. In the second region 8 < v < 13, there is only one kind of excitation spectrum which 
is unstable (i.e., u 2 < 0). In the third region v > 13, there exist three kinds of excitation 
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spectrum. One of them is stable, and other two are unstable. These unstable regimes are 
explicitly indicated in Fig. 2(b), which is the enlargement of the curves in the vicinity of the 
point A in Fig. 2(a). From Fig. 2(c), we see that u 2 has two positive real number solutions 
and one pure imaginary solution which is not indicated in the figure, these means that there 
does not exist solutions of oo 2 < 0. Therefore, the excitation spectrum of this case is always 
stable. 
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Figure Captions 

FIG.l. (a) Excitation spectrum uo 2 for the out-of-phase solution for 77 = 80, g — 7.29, 
and 7] m = 4.932. (b) The enlargement of the part of — 10 < uo 2 < 20 for the curve c. 

FIG. 2. (a) Excitation spectrum u 2 for the in-phase solution with = 4.584. (b) The 
enlargement of the curves in the vicinity of the point A in Fig. 1(a). (c) Excitation spectrum 
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uj 2 for the in-phase solution with £ m = 4.431. Here we set rj = 80, g = 7.29. 
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